Charge and Spin Ordering in the Mixed Valence Compound LuFe2C>4 
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Landau theory and symmetry considerations lead to a unified treatment of charge and spin order- 
ing in the mixed valence compound LuFe2C>4. The unusual evolution of charge ordering is attributed 
to interactions between charges and phonons whose frequencies are calculated from first principles. 

PACS numbers: 75.25.+z,75.10. Jm,75.40.Gb 



In the past few years there has been an explosion of 
interest in systems which display simultaneous magnetic 
and ferroelectric ordering. 1.2.3,4] In many of these com- 
pounds two magnetic phase transitions are required to 
produce the necessary lowering of symmetry to allow 
ferroelectricityQ. A slightly different scenario is pre- 
sented by members of the mixed valence family RFe204, 
where R is a rare earth, in which half the Fe ions carry 
charge +2e and half carry charge +3e and here we con- 
sider the case R = Lu. LuFe2C>4 (LFO) has a trigonal 
Rim crystal structure |6| and the Fe ions form triangu- 
lar lattice layers (TLL's) arranged in bilayers with two 
Fe ions per rhombohedral primitive unit cell, as shown 
in Fig. la. At temperatures above 500K, the valence 
electrons can thermally hop so that all sites appear to 
have charge 2.5e.§| As the temperature is reduced from 
500K, charge ordering (CO) correlations develop at wave 
vectors associated with "root 3" ordering within each 
TLL (see Fig. lb) and correspond to CO within each 
bilayer such that the dipole moment of all bilayers are 
parallel. (We call these "ferro" (F) configurations.) At 
T = Tqq = 350K [10illl| there is a continuous transition 
at which CO develops 3 dimensional (D) long range or- 
der without any long range magnetic order. [9[ However, 
it is striking that the long-range CO occurs in an "anti- 
ferro" (AF) configuration in which the dipole moments 
of adjacent bilayers are oppositely oriented. [§] Angst et 
at argue that the previous report of a spontaneous 
polarization[10] in this CO phase may be an artefact of 
the small applied electric field. As the temperature is fur- 
ther reduced, the spontaneous polarization P (in a small 
electric field) quickly saturates and is constant between 
320K and 240K[l(|. At T c = 240K there is a transi- 
tion in which a ferrimagnetic state appears which, at low 
temperature, has a moment per unit cell which is very 
nearly 1/3 of the moment that would occur if all spins 
were parallel |l2fll3|. At this transition P abruptly begins 
to further increase so that [l0| 



P(T) = P(T C ) + a(T c - T) 
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where (3' appears to be somewhat smaller than 1. At 
T = 175K there is a transition which may be of magnetic 
origin. [r| The major questions raised by these behaviors 
are (a) what does Landau theory say about the incom- 



mensurate order wave vector? (b) why are the dominant 
fluctuations for T > Tqo F whereas for T < Tqo the or- 
dered phase is AF? and (c) what form does the coupling 
between magnetic order and the polarization take and 
what restrictions does symmetry place on the direction 
of the polarization of a bilayer? 
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FIG. 1: (Color online) a) Fe ions in LFO with charge-charge 
interactions J„ indicated. The hexagonal (conventional) unit 
cell contains three bilayers configured so layers are stacked 
in the order ABCABC (A=red, B=black, C=blue) with two 
sites in the rhombohedral unit cell labeled "1" and "2" which 
are related by a center of inversion symmetry, b) root 3 
structure of charge or spin ordering in a single TLL, where 
L — a cos ip, M — a cos(t/j + 2n/3), and N = a cos(tp — 2n/3), 
where a is the amplitude and ip the phase of the order pa- 
rameter. The shaded rhombus is the broken symmetry unit 
cell. 



We first consider question (a). A major advance in 
understanding CO was due to Yamada et a?[l4[ who pre- 
sented calculations which could explain the unique zig- 
zag X-ray pattern of spots at T = 270K[14j at vectors 
(in hexagonal rlu's) [l5[ 

qi = (/i + 1/3 + <5'/3,/i+ 1/3 + (573,7-) , 

q 2 = (h - 2/3 - 2*73, h + 1/3 + 573, t - 2) , 

q 3 = (h +1/3 + 573, h- 2/3- 2573, r - 1) , (2) 

where r = 3(n+ 1/2). Here h and n are integers and 6' — 
0.0081. However, surprisingly, Yamada et al. claimed 
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that 3D order would not occur even in the presence of 
the interactions J 2 and J3 between adjacent bilayers. 

We start with a Landau analysis of CO using the lat- 
tice gas model 14] in which one introduces a variable at 
the nth site of the rhombohedral unit cell at R, Q n (R), 
which assumes the value +1 (—1) if the site is occupied 
by a Fe 3+ (Fe+ 2 ) ion. As shown in Fig. 1, one has 
interactions J\ and J4 within each TLL, and interac- 
tions J 2 and J 3 between nearest neighbors on adjacent 
bilayers. J\ is AF (positive). As argued in Ref. 1X41 . 
rather than use a long-ranged Coulomb interaction jig, 
we invoke a short-ranged interaction strongly screening 
by the large dielectric constant [I?) of the bulk. The 
Fourier transformed variables, Qn(q), satisfy Q n (R) = 
J2q Qn{q) exp(— iq • R). At quadratic order the Landau 
free energy is 



1 2 

F = oE E Fnm(q)S n (qrS m (q) , 



(3) 



q n,m=l 



where S n (q) = (Q n (q)), where here and below (...) indi- 
cates a thermal average, and F nm (q) — F mn (q)* . Under 
spatial inversion, X, ZSi(q) = S*2(q)*, which implies that 
Fn = F 22 . In a simple approximation 

F nm (q) = cTS nm + J(Q, rn; R : n ) exp(iq • R) , (4) 

R 

where we set ks — 1, </(0, m; R, n) is the interaction be- 
tween sites m in the rhombohedral unit cell at the origin 
and n in the rhombohedral unit cell at R, and c is a 
constant of order unity. In Cartesian coordinates 

Fn = cT + 2Ji[cos(aq x ) + 2cos(aq x /2)cos(V3aq y /2)] 
+ 2J i [cos(^/3aq y ) + 2cos(3aq x /2) cos(V3aq y /2)] , 
F21 = J 2 e- 2 ^/ 3 A( qx ,q y y+J 3 e- lc ^/ 3 A(q x ,q y ), (5) 

where 

A{q x ,q y ) = 2e iq y a ^ /6 cos{aq x /2) + e ~ laq ^ /3 . (6) 

To identify the critical wave vector(s) that minimizes the 
eigenvalue of F we set aq x = 4ir/3 + rj and aq y = p. The 
wave vector for r\ — p = (which we refer to as the X 
point) gives the root 3 structure (see Fig. lb) in each 
TLL and it is expected that the critical wave vector will 
have small rj and p. To lowest sensible order, we find that 

F n = [-3 + (3/4)<5 2 - (V3/8)6 3 cos(30)]Ji 

+ [6-(9/2)<5 2 ]J 4 + 0(<5 4 ) , (7) 

where 77 = S cos <fi and p = 5 sin <j). Also 

A = -(V3/2)[t/ + *p] = (-V3/2)Je i *. (8) 

This is an important result which shows that generically, 
the X point (which can be stable for a 2D triangular 
lattice) is not stable in the present case due to the off- 
diagonal terms in A which are linear in the wave vector. 



This is an example of the Lifshitz condition which states 
that one can only have a temperature independent com- 
mensurate wave vector in a continuous transition if the 
wave vector has sufficiently high symmetry. (l8j 
Thus the critical eigenvalue, A(q) is 



A(q) = 
V3<5 



Jf + J 3 2 + 2 J 2 J 3 cos^ + cq z /3) 



1/2 



(9) 



We regard A(q) as a function of 0, S and q z = 2<f) + cq z /3. 
Minimizing with respect to 4> yields 4> — 4>n = 2nir/3, for 
n = 0,1,2. For A(q) to be minimal, the square bracket 
in Eq. © should be maximal. For J2J3 negative, the 
minima occur for cos(20„ + c^z/3) = —1 in which case 
one finds the wave vectors of Eq. @ as was found in Ref. 
14. 15] For J-j-h > 0, the minima occur for cos(20„, + 
ck z /3) = +1, in which case one finds the wave vectors 
given by Eq. ([2]) but with r = in. Since fluctuation 
corrections to mean field theory will be smaller at higher 
temperature, we use the scattering at high temperatures 
(dominantly r is integral) to deduce that J 2 J3 is positive 
and leave the CO phase data (r is half-integral) to be 
explained as a fluctuation effect. 

Minimizing A(q) for J 2 J3 positive yields 



I Ji + M/Ji = V3S = 47n57V3 = 0.064 



(10) 



Also for Fn to be minimal at S — requires that J\ — 6 J4 
be positive. These relations indicate that the interaction 
J(r) decays with distance much more rapidly than do 
bare Coulomb interactions. So we neglect J4 in compari- 
son to J\. We set A(q) = to obtain the mean-field value 
Tqo ~ 3Ji. Since quasi-2D fluctuations can be severe we 
adopt the estimate J\ = 40 mcV = 500 K. 

We now discuss question (b) concerning the competi- 
tion between F fluctuations and AF ordering. Note that 
A(q) depends on q z only through the small interactions 
J2 and J3 with < J3 <§; J 2 - (The separation for J 3 (J 2 ) 
is 6. 30 A (3. 15 A) which leads us to posit that J3/J2 is of 
order 0.1.) Thus Eq. JTU1) indicates that J 2 /J\ = 0.06 
and hence that J 3 = O.IJ2 = 0.3 meV. The difference in 
T co for F and AF CO, AT, is, from Eq. ©, 

AT = T co {q z = 0) - T co {q z = Stt/c) = VSSJ 3 (11) 

from which we estimate that AT = 0.06 J3 = 0.25 K. Ac- 
cordingly, it makes sense to consider a theory in which 
F CO at q z = 0, with order parameter (OP) ctb, com- 
petes with AF CO at q z = 3tt/c, with OP a a- Our 
conclusion that AT is small is consistent with the sug- 
gestion of Ang st et al. that a spontaneous polarization 
is observed[10( because of the presence of a small applied 
electric field is enough to favor F stacking of bilayers over 
AF stacking. We point out that to obtain a polarization 
the small electric field has to also favor a commensurate 
state with <5 = in preference to an incommensurate 
state. We thus consider a free energy of the form[19] 
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va A a z B (12) 
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FIG. 2: a) Mean-field phase diagram [T^] near the bicritical 
point of Eq. (|12[) for v > 0. Since the dashed trajectory 
in the disordered (P) phase is closer to the F phase than to 
the AF phase, F fluctuations dominate AF fluctuations in 
the P phase, b) A single CO bilayer with a mirror plane M 
perpendicular to the wave vector which is along i. Squares 
are Fe 3+ 's and circles are Fe 2+ 's. 



We propose the trajectory shown in Fig. 2a to explain 
dominant F fluctuations in the disordered phase but with 
condensation into an AF ordered state. This trajectory 
implies a temperature dependent renormalization of the 
quadratic coefficients in Eq. (jT3J) which favors AF fluctu- 
ations over F fluctuations. Quartic terms for a single OP 
arise from implementation of the fixed length constraint 
and are the same for F (erg) and AF (a a) fluctuations. 
So these terms lead to a renormalization of Tqo > but do 
not explain the slope of the trajectory in Fig. 2a. Ac- 
cordingly we analyze a free energy of the form 

V = a\o*\Y + (1/2) X ?Y 2 , (13) 

where Y is a noncritical variable and a is either the OP 
a a or o~b- Minimizing with respect to Y we obtain a 



renormalized quartic interaction 

V = -(l/2) XY a 2 \a 2 \ 2 



(14) 



To leading order in the fluctuations we replace \a 2 \ 2 by 
4(|er 2 |)|(7 2 | and V leads to a correction to r of 



Sr = -4 XY a 2 (\a 2 \) 



(15) 



which is temperature dependent, increasing in magnitude 
as the temperature is lowered. It is crucial to show that 
such terms favor AF fluctuations and thus lead to a tra- 
jectory with the slope shown in Fig. 2a. 

Accordingly, we treat charge-phonon interactions 
where the phonon amplitude u is the noncritical vari- 
able and from Eq. (|13|) we see that u must be an even 
symmetry phonon. To see that this interaction favors the 
AF configuration look at Fig. 3 where we show the net 
force on a given TLL (called TLLo) from the neighbor- 
ing TLL's. One see that the forces from the neighboring 
TLL's tend to cancel (add up) if the adjacent TLL's are 
in an F (AF) configuration. Table I shows the results 
of a first-principles [20j calculation which gave optical 
phonons at 1 1 meV and, for the transverse g type needed 
for this mechanism, at 20 meV. 

To estimate the effect of this coupling we write 
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FIG. 3: (Color online) a) The in-plane components of forces 
on sites L, M, and N in TLLo (as in Fig. 1) in the F con- 
figuration. The solid line vector forces are from the charges 
(solid circles) in the TLL above TLLo and the dashed line 
vector forces are from the charges (dashed circles) in the TLL 
below TLLo. Lhe charges at A, B, and C are negative, pos- 
itive, and positive, respectively. The larger dots and heavier 
lined circles are charges of twice the magnitude of the smaller 
dots and lighter lined circles, b) The net force on the sites as- 
suming the separations between all planes are the same. For 
the AF configuration the dashed forces are reversed and the 
resulting total forces are listed as "Antiferro:". 



TABLE I: Calculated zone center phonon energies (in meV) 
and their symmetry labels. The F u and T g modes are odd 
and even under inversion, respectively. The nondegenerate 
(A) and doubly degenerate (E) modes correspond to motion 
along the c-axis and within the afe-plane, respectively. 



r 

Energy 


E U Eg A2u A\g A'2 U E u 

11.41 19.99 20.02 31.48 38.46 41.20 


r 

Energy 


Alg Eg Alu E U Eg A\g 

53.25 54.73 57.69 58.79 62.62 84.96 



where J(ry) is the heavily screened Coulomb interaction. 
As in Fig. 1, the Q's are given in terms of the order 
parameter ax, where X indicates either the F or AF 
configuration of TLL's. Because the transverse motion 
of planes is relatively soft, we consider displacements to 
lie within the TLLo plane. When minimized with respect 
to Ui, the free energy is 



Fph — 



z 2 \ax\ 4 
2Mu? D 



r) ( r ) ^ 



(17) 



where £ = [r/ J)(dJ/dr), rii is the component of r within 
the TLL, and z is the effective number of nearest neigh- 
bors. The actual number of nearest neighbors is 6, but 
since the forces do not all add up, we take z = 3 for the F 
configuration and z = for the AF configuration where 
the forces from adjacent TLL's nearly cancel. We set 
hujjj — 20 meV, r = 5 A, rii = 2 A, and use our previous 
estimate that J — 1 meV. The factor £ depends on how 
rapidly the interaction decreases with distance. For bare 
Coulomb interactions £ = — 1. But we are far from that 
regime. We take £ = —10, which is a value often found 
in insulators^i Also the Fe mass is M = 60 amu. Thus 



Fph 



-FoK F \ 



(18) 
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where F « O.lmeV. Equations (JT4J) and (fl5|> . yield 

<5r AF « -lmeV(|a 2 F |) , (19) 

which is enough to shift ordering from F to AF below T c . 

We now analyze magnetic ordering using a Landau free 
energy assuming only the z-components of spin are rele- 
vant. Symmetry dictates that the free energy quadratic 
in the magnetic variables, S n (q), has exactly the same 
form as Eq. ([3]). The important difference from CO is 
that for spin ordering one has significant reflections only 
at integer values of ck z /(2ir). Therefore in Eq. © it is 
clear that J2J3 is positive, so that the minima are given 
by Eq. @ with r = 3n. Although no such zig-zag pat- 
tern has been reported, such a phenomenon may be hard 
to detect if 8' is small, as it was overlooked by many 
studies of CO. We note that this spin incommensurabil- 
ity need not be coupled to that of CO. 

Now we consider the situation when the incommensu- 
ration 8 is zero. At low temperature CO is such that two 
thirds of the sites in a TLL unit cell are Fe 2+ and one 
third Fe 3+ , or vice versa. For spin ordering two thirds 
of the spins should be fully up and one third fully down. 
Note that the commensurate state with a single Fourier 
component is not satisfactory because it does not permit 
the expected nonzero total charge or magnetic moment 
in each TLL. To get that we invoke the interaction 

V = - a( 7(q) 3 a(0,0,3^) + (l/2) X - 1 a(0,0,3^) 2 (^0) 

where a is a constant. This term is allowed because 
3q + 3q z k is a reciprocal lattice vector. Here we assume 
that q z for the order parameter a is either F or AF, so 
that 3cq z /(2ir) is either an integer or a half-integer. Min- 
imization of V gives a(0, 0,3q z ) = axc(q) 3 . It is a l so 
worth noting that the phase tp of the OP (see Fig. lb), 
although not fixed by scattering data, is fixed by a sixth 
order term in the free energy, 

V 6 = C[o-(q) 6 + <r(q)* 6 ] = 2C| ( t| 6 cos(6t/>) , (21) 



where C is a constant and cr = \a\ exp(iip). For C < 0, 
Vq is minimal for ip = mr/3, where n is an integer. 

We now discuss question (c), starting with the exper- 
imental result of Eq. (1) for P(T) measured in a small 
electric field. flG] Since P(T) is already ordered, the mag- 
netoelectric free energy is of the form 

F = (l^ix^SPf +xl 1 8Pl]+a<j s 2 [P(a s )} 2 (22) 

where Pm = P(a s = 0), 0-5 is the ferrimagnetic OP, x 
is the dielectric susceptibility, a is a constant, and 8P = 
P(cr<j) — P(as = 0). This free energy is invariant under 
time reversal (as — * — as) and inversion (P — * — P). 
Minimizing with respect to <5P we find that 8P± = and 

P|| (a s ) = P(a s = 0) + 2a X \\P(a s = 0)a% , (23) 

so that if a s ~ (T- T c f , then 8P\\ ~ (T - T c ) 2 ^ . This 
result comports quite nicely with experiment, since this 
exponent at T c appears to be close to 1. 10] 

Finally, we comment on the direction of P for a single 
bilayer. As mentioned, the polarization is zero unless 
8 = (the system is commensurate). In Fig. 2b we show 
a bilayer with commensurate root 3 CO. Symmetry under 
the mirror operation perpendicular to the wave vector q 
forces P to be perpendicular to q. Although calculations 
at zero temperature [9] indicate that P is in the direction 
from a negative charge in one TLL to a positive charge 
in the other TLL, this special result need not hold for 
all temperatures and it would be interesting to actually 
measure the direction of P as a function of temperature. 

In summary: in this Letter we have provided a phe- 
nomenological explanation of the charge ordering, spin 
ordering, and ferroelectricity in LUFC2O4. 

We thank A. D. Christianson and M. Angst for helpful 
correspondance, A. Boothroyd for introducing us to this 
subject, and E. J. Mele and T. C. Lubensky for several 
stimulating discussions. 
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